ChinaXiv 合 作 期 刊 


Vol 14 No.1 天 文 研究 与 技术 第 14 卷 第 1 其 
Jan., 2017 ASTRONOMICAL RESEARCH AND TECHNOLOGY 2017 年 1 月 


CN 53-1189/P ISSN 1672-7673 


实测 中 激光 测 距 方 程 的 完善 


9 1,2 ”之 人生 1 1,2 Ztb J 刀 zl 
程 东 升 “ ， 汤 傅 峰 ， 李 春晓 ”， 能 焰 恒 
(1. 中 国 科 学 院 云南 天 文 台 ， 云 南 昆明 650011; 2. 中 国 科 学 院 大 学 ， 北 京 ”100049) 


摘要 : 考虑 望远镜 存在 跟踪 随机 拌 动 及 大 气 湛 流 引起 的 光斑 拌 动 等 因素 ， 推 导出 新 的 激 
光 测 距 方 程 。 采 用 数值 模拟 的 计算 方法 ， 研究 分 析 激 光 发 散 角 、 望 远 镜 跟 踪 随 机 拉动 和 大 气 
潢 流 引 起 的 光斑 持 动 对 回 波 强度 的 影响 ， 得 出 望远镜 跟踪 持 动 及 激光 发 散 角 是 影响 测 距 系统 
性 能 的 重要 因素 。 
关键 词 ， 望 远 镜 跟 踪 持 动 ， 大 气 沸 流 ; 激光 测 距 方程 ; 激光 发 散 角 
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激光 测 距 方 程 是 研究 激光 测 距 系 统 探测 能 力 的 基本 理论 依据 ， 对 激光 测 距 方 程 进行 分 析 能 够 对 测 
距 系 统 中 使 用 的 激光 融 的 能 量 、 激 光 发 散 角 等 参数 的 最 优选 取 提供 理论 支持 。 传 统 激光 测 距 方程 的 推 
导 过 程 中 假设 激光 脉冲 能 量 的 空间 分 布 是 均匀 的 ， 在 激光 大 气 传播 过 程 中 仪 考虑 了 大 气 的 衰减 作 
用 中 。 然 而 ， 在 激光 测 距 系统 中 ， 激 光 能 量 的 分 布 形式 ， 大 气 汕 流 对 激光 传输 时 抖动 、 扩 展 、 闪 烁 
以 及 望远镜 在 实测 时 的 跟踪 拌 动 都 会 导致 回 波 强度 的 变化 。 文 [2] 在 传统 激光 测 距 方 程 的 基础 上 ， 考 
TT 虑 了 激光 束 本 身 的 分 布 形 式 和 大 气 消 流 对 激光 传输 的 影响 ， 推 导出 激光 测 距 方程 的 新 的 表达 形式 ; 文 
[3] 基 于 对 望远镜 跟踪 抖动 满足 高 斯 分 布 的 假设 ,推导 出 含有 望远镜 跟踪 拌 动 参量 的 激光 测 距 方 程 。 
在 实测 过 程 中 ， 大 气 油 流 引起 光束 的 拌 动 、 扩 展 、 闪 烁 以 及 望远镜 的 跟踪 拌 动 同时 影响 回 波 强 度 ， 所 
以 需 将 大 气 渍 流 与 望远镜 跟踪 拌 动 同时 考虑 到 激光 测 距 方程 中 。 


1 激光 能 量 为 高 斯 (Gaussian) 分 布 时 能 量 表达 式 


在 激光 测 距 系统 中 使 用 的 激光 器 ， 甚 光 谐振 腔 发 出 的 基 模 辐射 场 的 横 截 面 振幅 与 能 量 是 高 斯 分 
布 。 激 光束 在 传输 距离 为 = 时 能 量 分 布 BCp) 和 截面 半径 |p|=p(z，y) = Ve + 入 之 间 的 关系 为 

二 (1) 
其 中 ,4 为 系数 ; w(z) 为 表征 激光 传输 距离 z 时 的 光斑 半径 


入 


En (2) 
TO 


wo (2) =w0 [1+( 


其 中 ，wo 为 激光 的 束 腰 半径 ;z 为 激光 传输 的 距离 。 
假设 光束 中 心 能 量 的 1/e 处 所 对 应 的 光束 半径 为 p,， 故 有 


3 


Pe 
= 3 
ww (z)” “39 
将 (3) 式 代入 (1) 式 : 2 
E(p)= Ae wo2. (4) 
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对 (4) 式 进行 积分 可 求 出 激光 束 能 量 分 布 天 (p ) 与 整个 激光 脉冲 能 量 的 关系 式 : 


E, = eCp)pdpaep = | eh = mAp,” 
0 0 


Eo 2 
那么 ， E(p)= 5e oa ， (5) 


由 (5) 式 可 知 ， 当 空间 碎片 位 于 高 斯 光束 中 心 时 反射 的 回 波 强度 最 强 ， 当 空间 碎片 相对 高 斯 光束 中 心 
存在 偏差 时 ， 回 波 强度 降低 ， 激 光 能 量 为 高 斯 分 布 时 的 示意 图 如 图 1。 


图 1 激光 能 量 为 高 斯 分 布 时 的 示意 图 


Fig.1 Spatial distribution of the laser energy: Gaussian type 


2 ”望远镜 跟踪 拌 动 导致 目标 偏离 光束 中 心 距离 的 概率 分 布 表 达 式 


假设 望远镜 的 跟踪 拌 动 $ 是 满足 高 斯 概率 分 布 的 随机 变量 ， 此 时 $ 可 表示 为 
一 一 1 

g($; >) = 一 -一 

a | De 

其 中 ,$=(4,, 中 ,) ; 四, 为 望远镜 跟踪 抖动 由 在 方位 轴 上 的 分 量 ; 由 为 望远镜 跟踪 抖动 $ 在 高 度 轴 

上 的 分 量 ; 人 = (上 ， 几 ) 为 望远镜 的 固定 指向 偏差 ; 人 为 望远镜 方位 轴 上 的 固定 指向 偏差 ; jw, 为 望 远 


镜 高 度 轴 上 的 固定 指向 偏差 。 
当 望 远 镜 不 存在 固定 指向 偏差 时 ， 有 


03 AE NG) (6) 


gl$s 0 了) = 一 一， (7) 
27 | > | 
其 中 , T 为 转 置 运算 符 ， 区, 为 4 协 方差 矩阵 ， 可 表示 为 
GS Ou 
人 2 (8) 


当 望 远 镜 的 跟踪 随机 误差 在 高 度 和 方位 轴 互 不 相关 时 有 
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5 . 0 
三 (9) 
2 
, 0 oo, 


其 中 ，o, 为 望远镜 在 方位 轴 的 跟踪 精度 ; e, 为 望远镜 在 高 度 轴 的 跟踪 精度 。 
假设 望远镜 跟踪 拌 动 为 高 斯 分 布 且 仅 考 虑 望远镜 跟踪 随机 误差 ， 由 望远镜 跟踪 抖动 导致 的 空间 日 
标 偏离 光束 中 心 的 距离 p' 也 是 满足 高 斯 分 布 的 随机 变量 , p' 与 $ 的 关系 可 表示 为 


p. = 9R, (10) 
R 是 望远镜 与 被 测 目 标的 距离 。 所 以 有 
中 =D MAR. (11) 
将 其 代入 (7) 式 可 推出 . 
gp 0 了 (12) 
"| 


， [oR 0 
其 中 , 之 ，= 2 p2 | ” 


0 oR 


3 ”大 气 注 流 作用 导致 的 目标 偏离 光束 中 心 距离 的 概率 分 布 表 达 式 


由 于 大 气 滑 流 的 影响 ， 激 光 在 传播 过 程 中 不 断 随 机 地 改变 其 波束 特性 ， 表 现 为 光束 的 扩展 、 光 斑 
漂移 及 脉冲 宽度 的 增加 等 40 。 光 斑 漂 移 导 致 空间 碎片 相对 于 光束 的 中 心 距 离 六 随机 变化 ,假设 随 机 
变量 人 满足 高 斯 概率 分 布 函数 & ，g 可 表示 成 : 


si 0, 3 ,)= 


1 
on | | 
其 中 , p =(x, y,); ww 为 p 在 望远镜 方位 轴 的 分 量 ; y, 为 p 在 望远镜 高 度 轴 的 分 量 ; T 为 转 置 运算 


符 ; 0 代表 光束 的 中 心 为 随机 变量 的 均值 ， 区 , 为 太 的 协 方差 矩阵 。 当 大 气 汕 流 引起 的 光斑 漂移 在 
望远镜 方位 和 高 度 轴 互 不 相关 时 有 
Oo 0 
Z.-| a 


其 中 ，o; 为 望远镜 方位 轴 的 光斑 漂移 方差 ，o ,为 望远镜 高 度 轴 的 光斑 漂移 方差 。 
根据 文 [4] ， 大 气 庙 流 导致 光束 中 心 漂移 的 均 方 差 为 


i =WW.=% 
-5(pa) TY, (pa) (13) 


2 
Py 10. 0 | (14) 
kros D3 
当 光 斑 漂移 方差 在 望远镜 方位 轴 和 高 度 轴 相 等 时 ， 可 推出 : 
2 
02 = (15) 
pr D3 


其 中 ,为 大 气相 干 长 度 ;= 了 ;入 为 激光 波长 。 


4 综合 考虑 两 随机 作用 源 时 的 激光 测 距 方程 推导 


漫 反射 目标 相对 于 光束 中 心 的 偏差 距离 是 望远镜 跟踪 抖动 与 大 气 消 流 共同 作用 的 结果 , 且 这 两 个 
作用 相互 独立 ， 根 据 概率 论 ， 两 个 物理 过 程 对 同一 参数 的 作用 可 表示 成 两 个 函数 的 卷 积 “%”， 即 
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人 (16) 
具体 形式 可 表示 为 
1 i 二 党 
记 Pas 0， 六 ee ) = 一 一 Pu) ?9 (17) 
27 | by total | 
其 中 ， 可 表示 为 
oR +o? 0 o2 0 
Xl1 X2 at 
= + 三 二 ， (18) 
p> total > 12 2 遇 0 oR 十 是 | 0 同 


其 中 po 为 综合 考虑 两 个 作用 元 时 ， 漫 反 射 目 标 相 对 于 光斑 中 心 的 距离 偏差 ; T 为 转 置 运算 符 。 
ep ; ， 能 够 求 出 在 两 个 随机 作用 源 共 k 同 作用 时 ， 激 光 在 传输 至 z 处 的 能 量 均值 5. 


E, 
根据 E(p) = 
Tp。 


J a 35) dvdy 


J 


， 激光 发 散 角 半角 9.， 温 反射 目标 距 测 站 距离 为 RR 在 距离 为 R 时 的 光 
斑 面 积 为 Sy， 透 过 率 7。， a h 为 普天 元 常量 ; v 为 光波 
ee 透 过 率 ; a 为 大 气 衰减 系数 。 被 测 目标 为 均匀 漫 
反射 体 (反射 光 强 度 在 2 立体 角 内 各 向 同性 ) ， 有 效 反 射 面 为 S$， 反 射 率 为 p,， 当 5S 和 < 5% 时 ,口径 为 
D 的 望远镜 接收 光电 子 数 n, 可 表示 为 


(19) 


a ， DD 
ESp,T.T Tnm 


n, = (20) 
, 2TR hv 
将 代入 可 得 5 
风 Sp,T .TT nm 
0 4 
(21) 
1 /fp 2 p 2 27R hv 
.|y 下 四 区 四 区 
tota Oo Oo, 
其 中 ,，p. 为 光束 传输 距离 为 z 时 的 光斑 半径 ， 可 表示 为 
D 
Pp. = 7 + OR+p,., (22) 


po 为 大 济 流 引起 的 光 来 扩展 的 均值 ,假设 大 忆 济 流 引 起 的 光 来 扩展 的 均值 为 0， 将 其 代入 帮 程 中 可 
得 空间 雄 片 漫 反 射 激光 测 距 方 程 : 


D? 
总 Sp,T.T IT nm 
(23) 


nn 二 


nV/D/2+OR) +205 MDI+OR) +T20 27Rhy 


5 数值 模拟 计算 分 析 


5.1 ”望远镜 跟踪 精度 对 不 同 激光 发 散 角 测 距 时 的 回 波 强度 影响 
图 2 显示 了 在 忽略 大 气 消 流 对 光束 的 作用 时 ， 望 远 镜 口径 为 1.05 m， 日 标 距 离 为 1500 km， 激光 
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发 散 角 为 1” 时 的 回 波 能 量 分 别 与 激光 发 散 角 为 5”"、10” 时 回 波 能 量 的 比值 随 望 远 镜 跟 踪 精 度 的 关系 。 
A eo a ed 对 提高 回 
波 强度 的 作用 非常 明显 。 而 随 着 望远镜 跟踪 性 能 的 下 降 ， 减 小 激光 发 散 角 并 不 能 提高 望远镜 接收 的 回 
波 强度 。 


25 
dl: 一 激光 发 散 角 1 "的 回 波 能 量 与 10” 回 波 能 量 的 比值 随 望 远 镜 跟 踪 精 度 的 关系 
\ 一 激光 发 散 角 1" 的 回 波 能 量 与 5" 回 波 能 量 的 比值 随 望远镜 跟踪 精度 的 关系 


回 波 强 度 比值 


望远镜 跟踪 精度 /” 
图 2 激光 发 散 角 为 1"、5”"、10” 时 的 回 波 能 量 比值 随 望 远 镜 跟 踪 精 度 的 关系 


Fig.2 Echo energy ratio variations as a function of tracking precision with a laser divergence angle of 1”, 5” and 10” 


在 不 同 大 气相 干 长 度 条 件 下 ,激光 发 散 角 为 1” 时 ， 望 远 镜 跟 踪 精 度 对 回 波 强 度 的 影响 。 图 3 显 
示 了 在 激光 发 散 角 为 1， 望远镜 跟踪 精度 为 1 与 跟踪 精度 为 5”、10” 所 接收 的 回 波 能 量 比 值 随 大 气相 
干 长 度 的 变化 关系 。 由 图 ， 在 不 同 跟踪 精度 下 的 回 波 能 量 比值 随 大 气相 干 长 度 的 增加 而 提高 ， 在 大 气 
相干 长 度 为 1cm 时， 跟踪 精度 为 1" 与 10" 所 接收 的 回 波 能 量 比值 约 为 30， 而 当 大 气相 干 长 度 为 10 em 
时 ，1” 跟踪 精度 的 回 波 强度 比 跟踪 精度 为 10” 时 提高 了 近 80 倍 。 


100 
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地 60 

Pe 

害 ( 望远镜 跟踪 精度 为 1” 与 5” 时 的 回 波 强 度 比 值 随 大 气相 于 长 度 的 变化 关系 
-- 望远镜 跟踪 精度 为 1” 与 10” 时 的 回 波 强度 比值 随 大 气相 干 长 度 的 变化 关系 
伪 40 


20 


0 0.05 0.1 0.15 02 0.25 03 0.35 0.4 0.45 0.5 
大 气相 干 长 度 /m 
图 3 望远镜 跟踪 精度 为 1"、5”"、10” 时 的 回 波 能 量 比值 随 大 气相 干 长 度 的 关系 


Fig. 3 Echo energy ratio variations as a function of atmospheric coherence length with tracking precision of 1”, 5” and 10” 


5.2 大 气 漠 流 强度 对 不 同 激光 发 散 角 测 距 时 的 回 波 强 度 影响 
图 4 显示 了 在 忽略 大 气 满 流 作用 的 条 件 下 ， 激 光 发 散 角 为 1” 时 的 回 波 强度 分 别 与 发 散 角 为 5”、 
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10” 时 的 回 波 强度 比值 随 大 气相 二 长度 的 变化 。 由 图 ， 在 大 气相 干 长 度 为 1 cm 时 ， 激光 发 散 角 为 1 
时 的 回 波 强度 与 发 散 角 为 10” 时 的 回 波 强度 的 比值 约 为 10。 随 着 大 气相 干 长 度 的 增加 ， 比 值 也 随 之 提 
高 ， 在 大 气相 干 长 度 为 10 em 时 ， 二 者 比值 约 为 55 ， 而 且 随 着 大 气相 干 长 度 的 增加 回 波 强度 还 有 一 定 


的 提升 空间 。 

而 大 气相 干 长 度 对 回 波 强度 的 影响 

70 

60 
划 50 
社 为 5” 和 10" 时 的 回 波 能 量 比值 随 大 气相 干 长 度 的 变化 关系 
起 40 外 为 1” 和 10" 时 的 回 波 能 量 比值 随 大 气相 干 长 度 的 变化 关系 
党 
可 30 

20 

10 

0 

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
大 气相 干 长 度 /m 


图 4 激光 发 散 角 为 


1"、$"、10" 时 的 回 波 能 量 比值 随 大 气相 干 长 度 的 关系 


Fig.4 Echo energy ratio variations as a function of atmospheric coherence 


length with a laser divergence angle of 1”, 5” and 10” 


在 不 同 望 远 镜 跟 踪 精 度 的 条 件 下 ， 激 光 发 散 角 为 1 时， 大 气相 干 长 度 对 回 波 强度 的 影响 。 图 5 


显示 了 激光 发 散 角 为 1" 时， 大 气相 干 长 度 为 20 cm 的 回 波 能 量 与 相干 长 度 分 别 为 1cem、5 cm 和 10 cm 


时 ， 回 波 强 度 的 比值 随 望 远 镜 跟 踪 精 度 的 关系 。 如 图 ， 在 望远镜 跟踪 精度 为 1” 时 ，20 cm 与 1 cm 的 
大 气相 干 长 度 所 对 应 的 回 波 强 度 比值 为 2.8， 随 着 望远镜 跟踪 精度 的 降低 ， 大 气 油 流 强度 的 改善 并 不 
能 有 效 提高 回 波 强 度 ， 在 同样 的 跟踪 精度 下 ， 大 气相 干 长 度 从 5 cm、10 cm 分 别提 高 到 20 cm 后 ， 回 


波 强 度 分 别提 高 了 1.2 和 1.1 倍 。 


8 激光 发 散 角 为 1”, 不 


同 大 气相 干 长 度 下 回 波 能 量 比值 随 望 远 镜 跟踪 精度 的 变化 关系 


度 为 20 cm 与 1 cm 回 波 强度 的 比值 随 望远镜 跟踪 精度 的 变化 关系 


度 为 20 cm 与 5 cm 回 波 强度 的 比值 随 望远镜 跟踪 精度 的 变化 关系 
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一 大 气相 干 长 
一 - 大 气相 于 长 

一 大 气相 干 长 


回 波 能 量 比值 


度 为 20 cm 与 10 cm 回 波 强度 的 比值 随 望远镜 跟踪 精度 的 变化 关系 


望远镜 跟踪 精度 /” 
图 5 大 气相 干 长 度 为 1 cm、5 cm、10 cm、20 cm 时 回 波 强度 的 比值 随 望 远 镜 跟踪 精度 的 关系 
Fig.5 Echo energy 


ratio variations as a function of tracking precision with 


atmospheric coherence length of lcem, Scm, l0cm and 20cm 
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6 总 结 
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激光 测 距 系统 的 回 波 强 度 受 激光 发 散 角 、 望 远 镜 跟踪 精度 与 大 气 潮流 强度 共同 作用 的 影响 。 对 低 
轨 的 目标 进行 测 距 时 ， 当 站 址 处 的 大 气相 干 长度 不 小 于 10 cm 时 ， 大 气 潮流 造成 的 光斑 拌 动 对 测 距 系 
统 探测 性 能 的 影响 不 明显 ; 在 保证 望远镜 跟踪 精度 的 前 提 下 ,减少 激 光 发 散 角 能 够 大 幅 提 高 测 距 系 统 
的 探测 性 能 ， 然 而 在 激光 器 光束 质量 确定 后 ， 减 少 激光 发 散 角 需 增 加 望远镜 口径 ， 提 高 大 口径 望远镜 
的 跟踪 精度 对 伺服 控制 系统 及 驱动 系统 提出 更 高 要 求 。 在 现 有 的 技术 条 件 下 ， 利 用 测 距 系 统 成 像 终 端 
对 目标 进行 闭环 跟踪 ， 通 过 实时 控制 高 速 压 电 摆 镜 以 修正 激光 发 射 方向 ， 采 用 这 种 方式 来 减 小 望远镜 
跟踪 精度 对 测 距 系 统 性 能 的 影响 。 


参考 文献 : 
[1]  ” 叶 叔 华 , 黄 焉 . 天 文 地 球 动力 学 [M]. 山东 :; 山东 科学 技术 出 版 社 , 2000. 


[2] Xiong Yaoheng, Feng Hesheng. Modification of laser ranging equation | R]// Proceedings of 
the 13th International Laser Ranging Workshop. 2002. 


[3] Degnan J J. Millimeter accuracy satellite laser ranging: a review [J]. Contributions of Space 
Geodesy to Geodynamics: Technology, 1993. 133-162. 
[4] Yura H T. Short-term average optical-beam spread in a turbulent medium [J|]. Journal of the 


Optical Society of America, 1973, 63(5).: 567-572. 

[5] ” 范 建 兴 , 杨 福 民 ， 陈 启 秀 . 卫星 激光 测 距 中 的 激光 脉冲 波形 变化 的 理论 分 析 及 数值 计算 
[J]. 中 国 科学 (A 辑 ), 2001, 31(1): 63-69. 

[6]  ” 施 雨 , 李 浪 武 . 概率 论 与 数理 统计 应 用 [M]. 西安 : 西安 交通 大 学 出 版 社 , 1997. 


An Improvement of Laser Ranging Equation in Practical Use 
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Abstract: In this article a new laser-radar ranging equation is deduced. This equation takes into 
consideration the tracking jitter of a telescope and the beam jitter caused by atmosphere turbulence. Numerical 
simulation is applied to the equation on echo intensity; this paper also explores in detail factors such as laser 
divergence angle, tracking jitter of telescopes and beam jitter caused by atmosphere turbulence. Analyses 
indicate that both tracking jitter and beam jitter are crucial factors for the performance of a ranging system. 
Key words: Tracking jitter of telescopes; Beam jitter caused by the atmosphere turbulence; Laser ranging 


equation; Laser divergence angle 


